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. n / Q 0 SUMMARY 

1 /b 

Pulsating aurorae are proposed as a source of the in- 
frasonic waves associated with geomagnetic activity re- 
ported by Chrzanowski et al. One of the most plausible 
generation mechanisms of these long period pressure 
waves is the periodic heating of the air around 100 km, 
corresponding to the auroral coruscation reported by 
Campbell and Rees. In order to show the energetic rela- 
tionship between source input and pressure change at sea 
level, some theoretical calculations are performed by 
using a simple model of auroral distribution in the iso- 
thermal atmosphere. _ , . 
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INFRASONIC WAVES FROM THE AURORAL ZONE 


by 

Kaichi X- Maeda 
Goddard Space Flight Center 

and 

Tomiya Watanabe 
University of British Columbia 


INTRODUCTION 

The purpose of this paper is to interpret the origin of the strange traveling atmospheric 
waves observed at the ground during intervals of high geomagnetic activity (Reference 1). In 
particular, discussion centers on possible modes of atmospheric oscillations caused by the 
periodic bombardment of auroral particles in the polar mesosphere. 

Trains of the waves are detected by a system of four microphones placed on each corner of 
a quadrant roughly 8 km square, located north of Washington, D. C. The presence of a traveling 
wave is established when the same wave forms can be found on all four records with certain 
time shifts between them. These time displacements are used to determine the direction of 
wave propagation and the horizontal phase velocity of the waves. The periods of these infrasonic 
waves are usually 20 to 80 sec, but occasionally 100 to 300 sec waves are recorded. The pres- 
sure amplitude ranges from about 1 to 10 dynes/cm 2 . 

One of the peculiar features of these waves is the change of arrival direction with time of 
day. The general trend is from the northeast in the evening, from the north about midnight, and 
from the northwest in the morning. The shift back to northeast is somewhat discontinuous. 

These waves occur less frequently in the daytime. 

Since auroral activity predominates at about midnight local time, the time dependence of the 
appearance of infrasonic waves during intervals of high geomagnetic activity can be explained by 
assuming that the source of this kind of wave is located somewhere in the auroral region, as can 
be seen from Figure 1. The left side of this figure shows the diurnal variation of the arrival of 
sound waves during magnetic storms, reported by Chrzanowski et al. (Reference 1). This as- 
sumption is consistent with the finding that the enhancement of the wave intensity at Washington, 
D. C. is generally delayed several hours after the increase of geomagnetic activity, except in the 
case of very severe magnetic storms. During these storms the delay is less than 30 minutes, 
indicating sourthward spread of the source of these waves near Washington, D. C. 
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In attempting to understand these peculiar pressure waves, it is proposed that they originate 
from a certain domain of the ionospheric region in the auroral zone which is heated periodically 
by a severe bombardment of auroral particles. It is thought that this periodic precipitation of 
auroral particles occurs simultaneously with geomagnetic pulsations. 

INTERRELATIONS BETWEEN GEOMAGNETIC FLUCTUATIONS, PULSATING 
AURORAE, AND INFRASONIC WAVES 

The auroral luminosity often fluctuates with incoming auroral particles. The fluctuations of 
these particles, mostly electrons, can be explained by either a periodic change in the acceleration 
mechanism of incident particles or by the change of mirror heights of trapped particles in the 
earth's atmosphere, following the variations of the field intensity of the earth's magnetosphere. 
The variations of geomagnetic field intensity and of auroral brightness are therefore closely re- 
lated. A clear example of the correspondence between the pulsating aurorae and the rapid varia- 
tion of geomagnetic intensity observed at the ground has been given by Campbell (Reference 2). 

According to the analysis of the space probe data given by Sonnet et al. (Reference 3), and 
that of Explorer X (1961 *1) given by Heppner et al. (References 4 and 5), the region near the 
geomagnetic equator beyond about 10 earth radii is occasionally greatly disturbed. Such dis- 
turbances may be propagated along the magnetic lines of force as hydromagnetic waves, and are 
transformed into electromagnetic waves when they reach the conducting ionosphere (References 
6-8). After penetrating the ionospheric region as electromagnetic waves, they are almost per- 
fectly reflected at the earth’s surface (Reference 9). This produces standing hydromagnetic 
waves along magnetic lines of force. The simplest mode of such standing waves is the funda- 
mental mode, whose unique node is on the geomagnetic equatorial plane, with two loops of oscil- 
lation on the ends of the line of force on the earth, one in each hemisphere. This mode of oscil- 
lation has been investigated theoretically by Dungey (Reference 10), who called it the normal 
mode of torsional oscillations of the magnetic field in the earth's cavity. The eigen period of 
these oscillations increases rapidly with the latitude of the magnetic line of force intersecting 
the earth's surface, varying from several tens of seconds in the subauroral region to about 10 
minutes in the polar region (References 11-13). 

In addition to the above mode, there is another mode of oscillation responsible for geomag- 
netic pulsations of shorter periods. As noted by Dessler (Reference 14), the velocity of Alfven 
waves decreases very rapidly with decreasing height in the exospheric region below around 2,000 
km. Therefore, because of the continuity of energy flow Alfven waves of certain periods coming 
from the outer exosphere can be expected to be intensified. In this mode of oscillation, the layer 
of maximum Alfven velocity, approximately between 1500 and 3000 km, becomes the node of 
oscillation, and the loop of oscillation is near the earth's surface. This mode of hydromagnetic 
oscillation corresponds to the geomagnetic pulsations with periods from about 1 to several 
seconds, which appear more frequently in nighttime than in daytime (Reference 15). 

In considering the above facts, it might be speculated that a possible mechanism for the pro- 
duction of infrasonic waves during auroral activity would be the penetration of Alfven waves, 
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including modified Alfven waves and retarded sound waves (Reference 13), through the ionosphere. 
However, as will be shown later, these contributions are very small compared with the pressure 
disturbances produced by the periodic heating of the lower ionosphere caused by auroral 
bombardments. 

As shown by Heppner (Reference 16) auroral activity predominates around midnight local 
time and the active region extends towards lower latitudes with increasing activity; pulsating 
aurorae appear at this phase of auroral activity. In other words, among several types of auroral 
displays pulsating aurorae appear with the largest disturbance in energy and they occur in fairly 
low latitudes. 

According to Campbell and Ress (Reference 17) the peak of pulsating aurorae is around 100 
km, the bottom is at 90 km, the effective thickness is of the order of 20 km, and the most frequent 
period is from 6 to 10 sec. From the direct measurement of auroral particles by means of rocket 
borne detectors, the energy flux of auroral particles, mostly electrons, is of the order of several 
tens of ergs/cm 2 -sec at weak aurora and increases by more than a factor of 50 at bright aurora 
(Reference 18). Thus, the energy flux in a strong pulsating aurora can be estimated to be of the 
order of 10 3 ergs/cm 2 -sec or more. This figure is consistent with the estimate given by 
Chamberlain (Reference 19), based on measurements of auroral luminosity. 

On the other hand, the energy flux of hydromagnetic waves deduced from the magnetic pulsa- 
tion data is less than 10 ergs/cm 2 -sec below 200 km. The energy flux of hydromagnetic waves 
increases with height. However, the contribution to pressure waves in the lower atmosphere 
decreases with the increasing height of the source, as will be shown later. 

Another evidence for the present idea is the very good correspondence between the appear- 
ance of pulsating aurorae and that of infrasonic waves (Figure 2). The occurrence of pulsating 
aurorae shown on the right of Figure 2 is taken from the visoplot of auroral activity, reported by 
IGY World Data Center A, on days which correspond to the events, reported by Chrzanowski et al. 
(Reference 1), shown on the left of this figure. Two occurrences of infrasonic waves in 1957 were 
reported, but since no visoplots of auroral activity are available before the IGY, they are omitted 
from this figure. The arrival direction of these waves might deviate from their true direction 
because of strong wind systems in high altitudes. Therefore, the appearances of pressure waves 
do not necessarily match the pulsating auroral events in detail. However, it is quite clear that 
there is a close correspondence between infrasonic waves and pulsating aurorae. 


MATHEMATICAL TREATMENT 

It is obvious that ray theory is not applicable in this problem, because the wavelength is of 
the order of an atmospheric depth or more for the infrasonic waves observed. 

The excitation and propagation of long period pressure waves in the atmosphere have been 
investigated by several workers, mainly in two fields of geophysics, i.e., meteorology and ion- 
ospheric physics. The former consists of the study of atmospheric oscillations (References 20-24), 
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AURORAL VISOPLOT (pulsating only ) 
UNIVERSAL TIME 



Figure 2— Diurnal variation of the arrival direction of inffasonic waves 
during specific magnetic storms detected at the National Bureau of 
Standards, Washington, D. C. (Chrzanowski et al., 1960), and corre- 
sponding auroral activities observed in the northern hemisphere. 

and the study of micro-barometric disturbances (References 25-36). The latter is mostly con- 
cerned with the investigation of ionospheric disturbances and has been discussed by Martyn (Ref- 
erence 37), Sen and White (Reference 38), White (References 39 and 40) and Hines (References 
41 and 42). 

Although the subjects considered in the noted references are quite different, the mathematical 
treatments are essentially the same. They are based on a differential equation of velocity di- 
vergence, which is derived from three fundamental equations, i.e., the equations of motion, con- 
tinuity, and of first law of thermodynamics. The theoretical aspects of the problems in geo- 
physics, including oceanography, are reviewed in References 43-46. Mathematical technique has 
been developed to solve these problems, introducing the so-called field variables instead of simple 
hydrodynamical variables. 

In the present calculations, however, the classical method based on the equations of velocity 
divergence will be used, because of the convenience in comparing the results with observed data. 
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Notations and Fundamental Equations 

Preliminary notation will now be given: 


c 


D/Dt 

f 

g(0, “ g) 

k 

l 

P. P. T 

Po ’ Po * T o 
p , p , T 

p T , p T , t t 

q' (x, z) 


R 




U(u, w) 


X, z 

7 


V 

k 

cr 

T 

x( x > z) 
q 


velocity of sound in the atmosphere in cm/sec, c 2 = ygH where H = RT/g is the 
scale height of the isothermal atmosphere, 

Eulerian derivative, d/dt + UV, 

resultant of external forces except gravity, dynes/gm, 
acceleration of gravity, g = 980 cm/sec 2 , 

horizontal wave number corresponding to \, in cm” 1 , i.e., 2n/k y 
vertical wavelength of a pressure wave in cm, 

small departures from static values of pressure, density, and temperature, 
functions of x, z, and t, in dynes/cm 2 , gm/cm 3 , and °K, 

static values of pressure, density, and temperature, functions of z only, 

static values of pressure, density, and temperature of air at sea level, 

total pressure, density, and temperature, i.e., p 0 + p, p 0 +p, and T 0 +T, 

rate of net accession of heat, erg/gm-sec; s(x, z) = (y- i)q where q = p T q' in 
ergs/cm 3 -sec, 

gas constant of air, B/M = 2.87 x 10 6 ergs/gm-°C, where B is the universal gas 
constant, 8.314 x 10 7 ergs/mol-°K, and M is the molecular weight of air, approxi- 
mately 28.97, 

velocity vector, where u is the horizontal (southward) and w is the upward com- 
ponent of air flow in cm/sec, 

horizontal (southward) and vectical (upward) coordinates, 

ratio of the specific heats of air, c p /c v ~ 1-4, where C p and C v are the specific 
heat of air at constant pressure and constant volume, respectively, 

entropy of air in ergs/gm-°K, 

horizontal wavelength of a pressure wave in cm, 

frequency of the pressure wave corresponding to r, i.e., 2n/r y 

period of pressure wave in sec, 

the divergence of velocity in sec" 1 , i.e., x = du/dx + dWdz , 

Coriolis vector in rad/sec. 


In Eulerian notation, the equation of motion is 

i 

DU 1 

Dt + p V P + S + axlJ ” f 


( 1 ) 
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and the equation of continuity is 


_D_ 

DT 




0 


(2) 


The variations of pressure due to thermal excitations can be derived from the first law of 
thermodynamics, which gives the change of entropy (References 43-45), 


S? = * • 


By making use of the perfect gas function for entropy (Reference 46) it can be shown that Equa- 
tion 3 is equivalent to 


Dp 

TFt - c 


2 


ft 


s(x, z, t ) 


(4) 


where 


s(x, z, t) - (y-l)yoq' (x, z, t) = (7 - 1) q(x, z, t ) . (5) 

Since, in the present problem, the period of oscillation is less than a few minutes, the Coriolis 
force due to the earth's rotation is negligible, and all other external forces can be assumed to be 
zero. The equation of motion for the present problem is then simply 

DU 1 _ _ 

w + To Vp * • (6) 


One-Dimensional Model 

For the one-dimensional case atmospheric motion has only the z component of velocity, w; 
therefore, in its linear approximation form, the equation of motion (Equation 6) is 

t9w _ ~<2p 

Po tz - ~jr - ps • (7) 

The equation of continuity in its linear approximation form is 

dp (9w Po 

■Jt + Po a - TT w = 0 ■ (8) 

And the linear approximation form of the entropy equation (Equation 4) is: 

dp , c/w 

Tt ~ ~Po c 71 + Po e* + (y~ 1 )«i ■ (9) 
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where q = q(z, t) is the periodically changing heat source, which can be assumed to be 


q(z. t) 



for z > 0 , 


for - z 0 < z < 0 , 


( 10 ) 


where 


q(z) 




-2 /h 


h 


(ID 


h is a constant and q 0 is the maximum rate of heat generation in an atmospheric column with unit 
cross section, in ergs/cm 2 -sec; z = 0 is taken as the height of the base of heating; and the value 
at the earth's surface is given by z = - z 0 . 

By eliminating p and p from the above equations, the following differential equation is obtained: 


d 2 w 1 aw 1 d 2 W 

7^--fraz c z at 2 


y ~ i <?q 
p oC 2 JT 


( 12 ) 


In the region z >0, where the heat source exists, 


1 ik ^ 

w " k 2 - k, e 


f 

Jo 


e' ikl? Q(^)d^ + C 


+ k - k e 

Kj K 2 


T 


e ' ik ^ Q(£)df + D 


(13) 


where k, and k 2 are given by 


= 1 /— — 

V l* 4H 2 


- _ I fl L_ 

VI 2 4H 2 


(14) 


and 


1 cr 

7 = 7 • 

n/'\ - y “ 1 dc J<£> 

QK> " p„e> * 

C and D are integration constants, to be determined by the boundary conditions. The time- 
dependent factor, e iat , is dropped in Equation 13 and will be always dropped hereafter. 
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For the region - z 0 £ z < 0, where no heat source exists, 


1 ik, z . 1 ik ,z 

* = e A + 


k — v e 2 B ; 
K 1 K 2 


(16) 


where A and B are constants. The condition that w = 0 at the earth's surface gives 


B = e 


i( k 2 -k 1 >* 0 A 


(17) 


The vertical velocity w should be continuous at z = 0, and therefore 


A - B = C - D 


(18) 


The pressure p also should be continuous at z = 0, so that the following condition is to be satisfied: 


, / dw\ 

9 ' c Ku +e 


(y-l)q(O) = Pt 


, ( dw\ 

'« c wL= 0 -< 


(19) 


where • - 0 and 


P C C°o)z = 0 


(20) 


the atmospheric density at z = 0. 
The above equations give 


k 1 C-k 2 D-k 1 A + k 2 B - 


k, - k 2 


P„ C 


(y-l)q(0) 


( 21 ) 


For z -oo, there should be no wave downward, thus from Equation 13 


■ -r 

*0 


C = - I e _ikl ^ Q(0 d£ ■ 


( 22 ) 


From Equations 18 and 19, 


y — 1 

A = C r q(0) . 

P c c 2 


(23) 
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To find the pressure variation at the earth's surface, we must obtain (dw/dz) i= , which is given 
by using Equations 16 and 17; 



“ ik l z 0 

e 


A . 


Then Equation 9 gives the pressure variation at the earth's surface, 


(24) 


— ^ 2 i z n 

E o ‘ a P s c e A ■ 


(25) 


where p s is the atmospheric density at the earth's surface. By substituting Equation 22 into Equa- 
tion 23, 


(y-l)Qoi/ h h 2 V 1/2 


P c c d 


t \i - n + yj 


(26) 


provided that the angular wave frequency <x is larger than a critical frequency^ , 


_ yg 
a A - TE 


(27) 


By using Equations 14 


o/ 2H 


(28) 


Thus the final expression of the amplitude of pressure variations at the ground, caused by the 
periodic disturbance in the upper atmosphere q(z, t), is 


p s 




h 

TT + 



(29) 


The values for p s /i 0 which may be derived from this equation are plotted vs. h in Figure 3 
for several values of angular frequency, a - . The full lines and dashed lines stand for the scale 
heights of 8 km and 6.8 km, respectively, for an isothermal atmosphere. The latter gives the 
proper ratio of ^p c yp^j 1/2 as compared with the observed atmospheric densities. From this 
figure, it can be seen that periodic heating corresponding to a flux of the order of 100ergs/cm 2 -sec 
produces pressure waves with amplitudes of the order of 1 dyne/cm 2 at the ground. 
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Two-Dimensional Model 


Equation of Velocity Divergence 

From Equation 6 the equation of motion for 
the two-dimensional case is 




The equation of continuity (Equation 2) can be 
written in a first order approximation as 


dp d PQ 


8 10 


Figure 3— Relative intensity of pressure waves on 
the ground, produced by periodic aurora! heating, 
vs. the scale height of the heat source, for the one- 
dimensional model with an isothermal atmosphere. 


where * - *<x, z, t) is the velocity divergence. By using Equation 32, Equation 4 can be written 


7t - Po g* - Po ° 2 x + s 


Under the assumption that the time variations of u, w , p, p 9 and s are proportional to a factor 
e icrt , the following relations between u, w, and p are obtained from Equations 30-33: 


d / 9 v 1 ds 

dx~ l c2 X” 8 w j “ ~p7~Sx ’ 


w = ci m + g 5F - - (35) 

iop = p Q gw - c 2 p Q y + s . (36) 

By eliminating u, w, and p from these equations, the following differential equation for the velocity 
divergence X (x, z) is obtained: 


V72 1 /dc 2 \ dx g dc 2 d x a 2 

X c 2 \ dz ~ r 7 dz ~ c 2 a 2 [ dz + (r_1) gj 5x 2 + ~ 


1 g /dc 2 \ d 2 s or 2 ( dc 2 \ ds o' 2 

^ > o ( z)c 2 - T7* N" * v *)~ 5 * - 71 y2 s 


dx 2 P 0 C 2 


The derivation of this equation is given in Appendix A. 
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The Diagnostic Diagram 

If there is no thermal excitation, the right side of Equation 37 is zero, and the solution of this 
homogenous differential equation corresponds to the free oscillation of the atmosphere on the 
nonrotating earth. 

We will now consider pressure waves traveling horizontally in this flat atmosphere. If u, w, 
p, and p are proportional to a factor the following is the equation for the vertical change 

of x( a > z ) : 


d2 x 1 /(Jc 


dz 2 


d X 

rsr “'ygj w 


- 1,2 


k 2 + 


gk 2 ( dc 2 

+ ■ 


Since c 2 = ygH(z) , this may be written 


d2 x; i d x 

^7r + -H (H'-i)ur + 




L ) 


= 0 


(38) 


(39) 


where H' = dH/dz . 

For simplicity, consider the case of an isothermal atmosphere, where H' = 0. Equation 39 
becomes 


d 2 V dy 

^--2N W +M 2 y = 0 , (40) 


where the constant N = 1/2H and 


M 2 



+ 


k 2 g 2 
c 2 <r 2 


(r- i~) 


(41) 


This differential equation has the solution: 


y(o-, z) = e Nl (Ae - ^ 1 +Be + ^) 


(42) 


where A and b are constants, 


M 2 = N 2 - M 2 > 0 


(43) 


for the noncellular solution, and 


IJ. - irj , 


v 2 = M 2 - N 2 > 0 , 


(44) 
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for the cellular solution. As was shown by Pekeris (Reference 27) for the noncellular solution in 
an isothermal atmosphere (Equation 42) the term with e^ z must vanish, i.e., B = 0. Otherwise 
the kinetic energy of noncellular waves, which is proportional to p 0 (z)*- 2 , would diverge. Further- 
more, because of the condition that the vertical component of velocity must vanish at the ground, 
we get only two possible types of free oscillation, specified by the following characteristic 
relations: 

Type I: k = ^ , v = u = c,andM = 8(2 ~ y) , (45) 

2c 2 

Type II: k = -g- . v = |u = ^v,and/ 2 =^ • (46) 

where V = cr/k and U = dcr/dk, the phase velocity and group velocity of noncellular horizontal waves. 
Type I corresponds to Lamb's wave (Reference 20), which exists at all frequencies; the waves of 
type II can be propagated only above a critical value a c , 

_ g _/y 

a c - C Vt (47) 

(the corresponding period is about 250 sec), and the wave is dispersive. 

It should be noted that, because of the decrease of atmospheric density, the amplitude of the 
Lamb wave (pressure variation) decreases with altitude by a factor exp(- gz/£ 2 ) , preventing the 
propagation of waves in any except horizontal directions (Reference 46). On the other hand, for 
the cellular solution rj stands for the wave number in the vertical direction and Equation 44 is 
equivalent to 


T}* 


(cr 2 - a A 2 ) + k 2 c 2 (a B 2 -a 2 ) 


(48) 


where 


- iZ 

A 2c 


_ g(y- 1) 1/2 

°B C 


(49) 


Brunt's frequency, cr B , is the frequency of the vertical oscillation of a free air parcel in the 
atmosphere, changing adiabatically. This expression was also derived by Vaisala as a stability 
parameter of the atmosphere (Reference 46). The frequency a A can be called the atmospheric 
resonance sound frequency (Reference 47). 


The curve for t> 2 = 0 shown in Figure 4 (a plot of cr vs. k) consists of two branches, A andB. 
The curve A starts from the cr-axis at o- = a A and is asymptotic to the line C, which corresponds 
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ANGULAR VELOCITY <7(sec“l) 



Figure 4— Diagnostic diagram for an isothermal 
atmosphere with scale height H = 8 km (T Q = 273 °K). 

Intensity of Pressure Waves at the Ground 


to the solution of the noncellular Lamb wave, 
cr = kc. The other curve B passes through the 
origin and is asymptotic to the horizontal line 
cr ~ o" B . Every point in the two regions where 
t) 2 > 0, one surrounded by curve A and the other 
surrounded by curve B, leads to an eigen solu- 
tion of Equation 39. In the third region (be- 
tween the other two), where rj 2 < 0, only the 
points on the line cr = kc and a = /kg (a > <x c ) 
lead to eigen solutions. Figure 4 is called the 
diagnostic diagram of the isothermal atmos- 
phere, and the waves corresponding to the two 
domains in which rj 2 > 0 are named Sonic (mode 
A) and Thermobaric (mode B), respectively 
(Reference 46). Sonic and thermobaric waves 
are called acoustic and internal gravity waves, 
respectively, by Hines (Reference 42). 


Since the source of excitation s(x, z) is limited to a certain area inside the auroral zone, the 
right side of Equation 37 is not uniform with respect to x. Therefore, this differential equation is 
not in general separable with respect to the variables x and z, and the following assumptions are 
made in solving the equation: 


1. The atmosphere is isothermal, with scale height H, in km. 

2. The distribution of the heat source is uniform along the y direction (this is assumed at the 
beginning to reduce the problem to two dimensions), but it is limited horizontally in the 

x direction within ± (i.e., - \ 0 < x < \ 0 ) and is extended vertically above a certain height z 

(i.e., z > z 0 ). 

3. The time variation of the heat source is periodic with an angular frequency o-, and the same 
phase within the domain indicated in item 2 above. 


The first assumption reduces Equation 37 to the following form: 


g B 2 \ 2 X d 2 X l dx 
a 2 ) dx 2 + dz 2 H dz 



1 

Po ( z )c 2 



g \ d 2 S 
Cr 2 Hj dx 2 



1 ds 
+ H Hz 


Assumptions 2 and 3 can be written as 


s(x,z,t) - s(z) e icrt e[z -z 0 ) [<9(x + x 0 ) -^(x-^ 0 )] 


(50) 


(51) 
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where 


s(z) 


ll ~(*“ z o)/ h 
h e 


(r- l)q 0 

h 


exp 


~( 2 ~ z o) 

h 


and #(<f) is a unit step function of £, 


*{£) 


1 forf>0 , 

0 for £ < 0 


(52) 


(53) 


It should be noted that q 0 is the maximum rate of heat release in the atmosphere in an air column 
of unit cross-section (in ergs/cm 2 -sec). 

Since the atmosphere is assumed to be isothermal, the density of air in equilibrium at height 

z is 


P 0 (z) - P, , (54) 

where p m is the atmospheric density at the earth’s surface and H is the scale height. 

In order to solve Equation 50 under the conditions listed, the following Fourier transforms 
are applied with respect to x. 


X(k, z ) 


S(k, z) 



e lk * X( x ’ z ) dx 


e* 11 * s(x , z ) dx 


(55) 


(56) 


Since the heat source s(x, z) vanishes outside of the auroral zone, both s(x, z) and *(x, z) must 
vanish at x = ±00. Thus 


1 

t/2tt 


1 




iis 

dx 2 




dx 


dx 


- k 2 X(k, z ) 


-k 2 S(k, z) 


(57) 


(58) 
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With these transforms, Equation 50 is written 


d 2 X dx 

— - 2NdF + M 2 X = F(«,k) , 


( 59 ) 


where 


F(z , k) 


' Po ( z > 


’d 1 S 1 dS / 


zr*- Td7- ki ( 

1- 


(60) 


and N and M 2 are the same as in Equations 40 and 41. The solution of Equation 59 can be written: 


X(k, z) - e 


' r 

'I 2m J 


F( z ' , k ) e " l2 dz ' 


+ e 


C 2 + 




k) e 


dz 1 


(61) 


where F(z ' , k) is given by Equation 60 and n, and n 2 are the roots of the characteristic equation 


n 2 - 2Nn + M 2 = 0 


(62) 


We assume 


nj - N - m 

rij - N + m 


(63) 

(64) 


where 


M = Vn 2 - M 2 > 0 for N 2 > M 2 , 

M = iTj, 7 7 2 = M 2 - N 2 > 0 for N 2 < M 2 

The integration constants C x and c 2 are determined by the following two boundary conditions: 
1. The vertical component of velocity vanishes at the ground, 

w(x, z = 0) = 0 


(65) 


( 66 ) 


2. The kinetic energy of the waves at z = « is either zero or remains finite. In the latter case, 
the vertical component of the disturbance should never be propagated downward at z = 00 . 
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The pressure change at the ground can be obtained from X (x, z = 0), which is given by the 
inverse Fourier transform of X(k, z = 0), i.e., 


x(*. 0) 


l 



X(k, 0) dk 


( 67 ) 


where X(k, 0) is the solution of Equation 59 at z = 0 satisfying the above conditions. It is given 
by Equation B15 in Appendix B. 

Since w(x, 0)= 0, then s(x, 0) = 0 and Equation 36 reduces to 


P(x, 0) 


C 2 P, 

x(*. 0) 


( 68 ) 


where *(x , 0) is given by the inverse Fourier transform of X(k. 0). The evaluations of p(x, 0) are 
given in Appendix C and the results are shown in Figures 5-10. Various values are given in Tables 
1 and 2. The symbols p r and p ; signify the real and imaginary parts of the cellular mode, respectively. 



Figure 5 — A plot of the intensity of the cellular 
waves p r and p. on the ground, as a function of the 
input energy flux q 0 for r = 30 sec and h = H, vs. 
x, the horizontal distance from the center of the 
source. Full lines and dashed lines correspond to 
= 10 km and X. Q = 1 km, respectively. 



Figure 6 — The same plot as Figure 5 except that 
h = 1/2H instead of h = H. 
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H and h = 1/2H, respectively. 






ATTENUATION 


Because of the viscosity and thermal 
conductivity of air, acoustic waves in the 
atmosphere attenuate. The attenuation coef- 
ficient a (r), in cm' 1 , for a wave of period r 
is given approximately by (Reference 48) 


a(r) 


4tt 2 
r 2 


1 



7 “ 1 

y 



(69) 


where the coefficient of kinematic viscosity 
v is approximately 


1.7 x 10~ 4 


(70) 


in cm 2 /sec, and the coefficient of thermal 
conductivity a 2 is 


a 


2 


in cm 2 /sec. 


2.1 x 10' 5 

P( z ) 


(71) 



Figure 9— The same plot as Figure 7 except that 
r = 100 sec. 


Since the air density p decreases exponentially with height, a 2 and v increase exponentially 
with altitude. The so-called attenuation factor, 


1 a ( T ‘ z ) " ex P 


f 

Jo 


a (r, z ' ) dz * 


(72) 


is shown in Figure 11 as a function of height, for r = 10 sec, 30 sec, and 100 sec. 

It should be noted that the relative amplitude of the pressure wave grows as it propagates 
upward by a factor e l/2H , where H is the scale height, whereas the absolute amplitude decreases 
by a factor e* z/2H , because of the exponential decrease of air density (Reference 49). Therefore, 
if the amount of excitation energy is the same, the absolute intensity of the pressure wave at the 
ground increases with the base of the excitation level. This is shown by Equation 29 and Equa- 
tions C9 and CIO with a factor e z °^ 2H for p s , p c (z = 0) , and p nc (z = 0), and by a curved dash 
line in Figure 11. 

However, because of the steep increase of the kinematic viscosity of air with altitude, the 
yield of pressure wave excitation drops sharply above a certain altitude for a given period (or 
frequency) of the wave. This is shown in Figure 12 for wave periods of r = 10, 30, and 100 sec. 
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Figure 1 0— Intensity distribution of nonceilular waves on the ground, p nc , in the units of q Q/ vs. horizontal 
distance x from the center of the source of width k Q = 10 km. Notice that the vertical scales are magnified by 
a factor attached to each curve. Full lines and dashed lines stand for h = H and h= l/2H, respectively. 


Table 1 

Intensity of the Cellular Mode of Infrasonic Waves on the Ground Below the 
Center of the Source in Units of q Q . 


h = H 

Pr 

Pi 

Pc 

h = 1/2H 

Pr 

Pi 

Pc 

k Q ~ 1 km 

r = 10 sec. 
30 sec. 
100 sec. 

6.74 x 10' 6 
6.15 x 10 -5 
5.46 x 10' 4 

4.40 x 10' 7 
1.21 x 10' 5 
4.26 x 10' 4 

6.75 x 10" 6 
6.28 x 10' s 
6.92 x 10" 4 

k Q = 1 km 

r = 10 sec. 
30 sec. 
100 sec. 

1.33 x 10' 5 
1.10 x 10‘ 4 
2.77 x 10' 4 

1.79 x 10" 6 

4.79 x 10‘S 
1.02 x 10" 3 

1.34 x 10-5 
1.23 x 10 -4 
1.06 x 10* 3 

k Q = 10 km 

r = 10 sec. 
30 sec. 
100 sec. 

1.82 x 10' 5 
4.29 x 10' 4 
5.25 x 10-3 

1.20 x 10' 6 
8.57 x 10' 5 
4.11 x 10' 3 

1.83 x 10* 5 
4.38 x 10* 4 
6.67 x 10*3 

= 10 km 

r = 10 sec. 
30 sec. 
100 sec. 

3.62 x 10' 5 
7.91 x 10' 4 
2.68 x 10'3 

4.80 x 10“ 6 
3.32 x 10" 4 
9.85 x 10-3 

3.65 x 10-5 
8.58 x 10" 4 
1.02 x 10* 2 

= 100 km 

r = 10 sec. 
30 sec. 
100 sec. 

1.82 x 10' s 
4.84 x 10' 4 
1.43 x 10' 2 

1.20 x 10" 6 
1.01 x 10" 4 
1.14 x 10' 2 

1.83 x 10" 5 
4.94 x 10' 4 
1.83 x 10" 2 

k Q = 100 km 

r = 10 sec. 
30 sec. 
100 sec. 

3.62 x 10' 5 
8.92 x 10' 4 

7.62 x 10" 3 

4.80 x 10" 6 
3.73 x 10‘ 4 
2.72 x 10* 2 

3.65 x 10" 5 
9.68 x 10' 4 
2.83 x 10* 2 
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Table 2 


Values for p nc , P nc/ /|p c |, and x Q1 the Horizontal Distance of the first Nodal Lines from the Center 
of the Source (horizontal width assumed to be \ 0 = 10 km). 


h = H 

Pnc 

Pnc/ Pel 

X 0 (km) 

h = 1/2H 

Pnc 

Pnc/|Pc 

x 0 (km) 

r = 10 sec 

1.28 x 1<T 9 

7.0 x 10' 5 

0.85 

r = 10 sec 

1.11 x 10 ' 9 

3.03 x 10' 5 

0.85 

30 sec 

3.24 x 10~ 8 

7.4 x 10 ' s 

2.56 

30 sec 

2.73 x 10 ' 8 

3.18 x 10' s 

2.56 

100 sec 

1.10 x 10' 5 

1.66 x 10' 3 

8.31 

100 sec 

1.02 x 10 5 

9.8 x 10' 4 

8.30 


It can be seen that there is an effective height of excitation of atmospheric acoustic waves for a 
given period and that this height increases with the period. 


As mentioned in the introduction, the excitation of pressure waves by hydromagnetic waves is 
only effective above the height of auroral activities. Figure 12 shows that waves produced at such 

a height are greatly attenuated before they 

i a2 

^ reach the ground. Therefore, the contribution 

/ r= ioo sec of hydromagnetic waves to the excitation of 




Figure 1 1 — Attenuation factor f a (r, z) vs. altitude of the z 

source z in km and -In f o ( r , z) vs. z: 

Figure 12— Relative yield of the source to the infrasonic 
r z waves in the isothermal atmosphere 

-lnf 0 (T # z) = J a(r, z ') dz ' , y— - f r z 

0 % ex p J dz< 

a dashed line is also used to indicate the exponential r 0 L_ J o 

amplification factor In (p z p 0 ) 1/2 ^z/2H / where the vs. the altitude of the source z # where the scale height 

scale height H is assumed to be 8 km. H is assumed to be 8 km. 








internal atmospheric waves is practically negligible. A more quantitative discussion will be 
given in the next section. 

CONDITIONS FOR WAVE FORMATION 

According to previous calculations, a pressure variation at the ground of the order of 1 dyne/ 
cm 2 can be expected if the maximum rate of heat generation is of the order of 100 ergs/cm 2 -sec 
and the layer of periodic heating is around 100 km altitude, with a thickness less than 10 km. 
Electrons with an energy of the order of 100 kev will lose most of their energy within a layer of 
the order of 10 km thickness around the height of 100 km (Reference 19, p. 290). According to 
Chamberlain, the rate of heat generation is the order of 60 ergs/cm 2 -sec in a bright aurora. 
Therefore, if the rate of heat generation changes periodically with this order of amplitude, bar- 
ometric oscillations of the order of 1 dyne/cm 2 at sea level can be expected from the sources in 
the upper atmosphere. 

However, several other conditions must be satisfied in order that the energy brought into the 
upper air by auroral electrons can be converted efficiently to pressure waves in the atmosphere. 
At first, the time r l9 in which an electron arrives at auroral height from outer space, should be 
smaller than the period, r, of the waves concerned. If r x »r, the phase of the time variation of 
the source differs from place to place, and the resultant pressure wave originating from those 
different sources is weakened by superposition. 

Since the auroral electrons, whose velocity is of the order of 10 9 cm/sec, lose their energy 
within a layer of 10 km thickness, r t < approx, 10~ 4 sec. This is much smaller than the period of 
the acoustic waves considered here. 

As discussed by Hanson and Johnson (Reference 50), electrons impinging into the upper at- 
mosphere lose their energy mostly by inelastic collisions with neutral air particles until ~ 2 ev, 
which is the lowest excitation energy of atomic oxygen (*D state). The time r 2 for l D excitation 
collisions with oxygen atoms is of the order of 10' 3 sec at 100 km and of the order of 1 sec at 
400 km. 

Below 2 ev, the electrons lose energy in the upper atmosphere mainly by elastic collisions 
with ambient electrons. The time for those low energy electrons to equilibrate with ambient 
electrons, r 3 , is obtained from the expression for the rate of energy loss of a fast electron im- 
mersed in a thermalized plasma (Reference 50), 


„ 5.7 x 10 3 E 3/2 

r 3 * r 


(73) 


where E is the electron energy, i.e., E £ 2 ev, and N e is the electron concentration in cm' 3 . The 
equilibration time, r 3 , given by this expression is about 10' 1 sec at 100 km and of the order of 
10~ 2 sec above 350 km. Thus r 2 + r 3 is much smaller than the periods of the acoustic waves 
under discussion. 
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It should be noted that elastic collisions with neutral particles dominate over those with 
ambient electrons below the F2 maximum (around 300 km). The time constant r 4 for these elastic 
collisions consists of r(0) and r n 2 , the time constants for the loss of excess electron energy to 
atomic oxygen and to molecular nitrogen, respectively (Reference 50), and is given by 


T = 2.08 x 10 n E~ 1/2 

4 n(0) + 2.36n(N 2 ) ’ V 

where n(0) and n(N 2 ) are the concentration per cm 3 of atomic oxygen and of molecular nitrogen, 
respectively. For E = 2 ev, r 4 is of the order of 0.1 sec at 100 km and increases with height. It 
is of the order of 100 sec at 350 km. 

As shown in the previous section, the excitation of atmospheric pressure waves above 200 km 
is not important; the time constant r 4 does not destroy the condition of wave formation. In other 
words, below 200 km r x +r 2 +r 3 +r 4 < approx, r. 

Another condition necessary for wave formation is that the time constant for cooling in a 
certain domain of auroral activity must be much longer than the period of oscillation. If the initial 
temperature T 0 is assumed to be horizontally uniform within the domain of the source, -X 0 <x <k 0 , 
then the temperature at the center after t(sec) is 


where 



(75) 


<D(y) 



(76) 


and a 2 is the thermal diffusivity (coefficient of thermal conduction) of air. T a is the temperature 
outside the source. 


The time, t a , necessary to reduce the initial temperature difference between the inside and 
outside of the source region can be estimated from 


2 ]£*" 


(77) 


K 2 


Since the thermal diffusivity of air a 2 at a height of 100 km is of the order of 10 6 cm 2 /sec 
(see Equation 71), the minimum source width necessary to satisfy the condition for thermal 
oscillation is of the order of 30m for a period of 10 sec, and of the order of 100m for a period of 
100 sec. The source width considered in the present calculation is significantly larger than 
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these widths. In other words, the time constant for source cooling is sufficiently long for pres- 
sure wave production. 

OTHER POSSIBLE MECHANISMS 

Now the possibility of other suitable mechanisms for the generation of microbarometric 
oscillations during magnetic disturbances will be considered. 


Periodic Heating Due to the Absorption of Hydromagnetic Waves 

As stated in the introduction, hydromagnetic waves from the magnetosphere lose part of their 
energy in the lower exosphere and in the ionosphere. The remaining energy, which leaks through 
the ionosphere, is observed as fluctuations of the geomagnetic field intensity or geomagnetic 
pulsations. 

The rate of energy dissipation of hydromagnetic waves in the ionosphere increases with in- 
creasing frequency (References 6, 51, and 52). According to Watanabe, the dissipation of the 
hydromagnetic wave energy in the ionosphere is negligible for waves with periods longer than 
20 sec. In the auroral region, the amplitude of geomagnetic pulsations sometimes exceeds several 
tens of gammas (i.e., giant pulsations) and the periods are sometimes longer than several tens of 
seconds. The intensity of the incident wave is usually smaller for higher frequencies (Reference 
15). As shown by the power density of the small scale fluctuations of magnetic field intensity 
observed in the earth's magnetosphere between 5 and 15 earth radii (Reference 3), the frequency 
spectrum of hydromagnetic waves, which are regarded as the origins of geomagnetic pulsations 
observed at the earth's surface, is a decreasing function of frequency above the ionosphere. In 
any case, the upper limit of the incident wave amplitude may be taken as 100 y. The energy flux 
associated with these hydromagnetic waves is then of the order of several ergs/cm 2 -sec 1 , under 
the assumption that the Alfven wave velocity above the ionosphere is of the order of 10 8 cm/sec. 
The rate of heat generation by absorption of these hydromagnetic waves in the upper ionosphere 
has been estimated by several authors and their results are shown in Table 3. 


Table 3 

Rate of Heat Generation Due to Attenuation of Hydromagnetic Waves with a 100 y Incident Amplitude. 


Source 

Period of 
HM waves 
(sec) 

Rate of Heat 
Generation in an 
Air Column of 
unit cross section 
(erg/cm 2 -sec) 

Thickness of 
Heating Layer 
(km) 

Altitude of the 
Center of the Heat 
Generating Layer 
(km) 

Dessler 
(Reference 53) 

1 

1.3 

100 

170 

Francis and 

6.3 

0.7 

50 


Karplus 
(Reference 51) 

0.3 

1.6 

100 

■ 

Akasofu 

1 

9.7 

100 

225 

(Reference 52) 


0.17 

200 

240 


100 

0.028 

300 

250 
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The rate of heat generation is generally less than 1 erg/cm 2 -sec, except when the period is 
1 sec. Therefore the amount of heat generated by hydromagnetic waves penetrating the ionos- 
phere is smaller than that due to the auroral particles. Table 3 also contains estimates of the 
thickness of the layer of heat generation. In every case the thickness is much larger for hydro- 
magnetic waves than for auroral particles. 

As shown in Figure 12, the effective height of pressure wave generation by periodic heating 
of the upper air is limited, and it is lower for the shorter period. Therefore, generation of 
acoustic waves by the attenuation of hydromagnetic waves is less efficient than generation by the 
periodic heating due to auroral particles, even if the energy flux of the incident wave is increased 
to be the same as the latter. 


Pressure Waves Due to the Impacts of Auroral Particles 

A particle coming into the atmosphere loses its energy by transferring its downward momen- 
tum to the air particles. A pressure wave can be generated if the flux of the particle changes 
periodically with time. 

The upper limit of the pressure intensity, from this process, can be estimated by assuming 
that all the incoming particles would stop in a very short time within a very thin layer. By as- 
suming the average energy of incident electrons is 6 kev and the maximum flux is of the order of 
10 10 cm' 2 -sec" 1 , the maximum pressure exerted upon the thin layer of virtual shock absorber is 
estimated to be of the order of 4.10' 8 dynes/cm 2 . The corresponding intensity at the earth’s 
surface is of the order of 10“ 5 dyne/cm 2 , if all electrons stop at an altitude of 100 km. 

Penetration of Hydromagnetic Waves Through the Ionosphere 

As mentioned above, most geomagnetic pulsations are due to hydromagnetic oscillations in 
the exosphere. These oscillations are related to electromagnetic oscillations in the space be- 
tween the earth’s surface and the lower boundary of the ionosphere. Any oscillation mode in which 
the compression of atmospheric matter is involved gives rise to a variation in the density, and 
consequently a pressure variation. Therefore, it may be possible that geomagnetic pulsations 
and microbarometric oscillations have the same origin, i.e., hydromagnetic oscillations of the 
earth’s exosphere. 

A simple situation can be found at the geomagnetic equator, where geomagnetic pulsations 
appear sometimes in the N-S component of the geomagnetic field, with amplitudes of the order of 
several tens of gammas and periods of several minutes (Reference 54). The origins of this kind 
of pulsation presumably exist in the outer boundary of earth’s magnetosphere near the geomag- 
netic equatorial plane, and it may be propagated as a modified Alfven wave. A modified Alfven 
wave is a transverse wave with respect to changes in the electric and magnetic fields. On the other 
hand, it is a longitudinal wave when viewed as fluid motion (Reference 55). An electromagnetic 
wave as well as a pressure wave should be in the space between the earth’s surface and the lower 
boundary of the ionosphere. The energy flux associated with an incident modified Alfven wave is 
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roughly (b 2 /8 itt) V a> where V A is the group velocity of the modified Alfven waves and can be taken 
as the Alfven wave velocity at a higher portion of the ionosphere. B is the amplitude of the wave, 

i.e., the intensity of magnetic fluctuation above the ionosphere. 

By assuming conservation of energy flux, the upper limit of the amplitude of the pressure 
wave at the earth's surface, p e , can be estimated as follows: 


where c s is the sound wave velocity at the earth's surface, and p s is the air density at the earth's 
surface. By taking V A = 3.10 7 cm/sec (corresponding to daytime at the sunspot maximum 
activity) B = 30y, c s = 3.10 4 cm/sec, and p s = 1.25 x 10“ 3 gm/cm 3 , we find p e ~ 3 dynes/cm 2 . 
In this estimate the conversion factor between the energy of incident Alfven waves and that of 
secondary pressure waves is assumed to be unity. This factor must actually be very small be- 
cause of reflections and energy dissipations of incident waves at the upper part of the ionosphere. 
If it is not small, infrasonic waves could appear in the equatorial region during strong magnetic 
disturbances. Since the occurrence of aurorae in these regions is negligible, this might provide 
a method of direct detection of modified Alfven waves coming into the earth's atmosphere from 
the magnetosphere. 


CONCLUSIONS 

We have shown that one of the most plausible mechanisms for pressure wave generation 
during geomagnetic disturbances is the periodic heating of the polar ionosphere by auroral par- 
ticles, observed as pulsating aurorae. As emphasized by Campbell (Reference 56) the main energy 
source for this type of auroral activity is not necessarily incident auroral particles, but a flow of 
secondary electrons called the electro-jet. In this respect periodic heating by these intermittent 
electric currents is essentially the same as the so-called Joule heating discussed by Cole (Ref- 
erences 57 and 58). 

The following conclusions can be drawn from the present calculation: 

1. From Figures 3, 5-9, and 13 it can be seen that an incident energy flux of more than 100 
ergs/cm 2 -sec will produce acoustic waves observable at the ground, provided the periods 
are longer than about 10 seconds. 

2. The relative intensity of the pressure wave at the ground is higher when the heating is con- 
centrated within a thin layer than when it is distributed over a wide range of altitudes 
(Figures 3, 5, and 6). 

3. The intensity ratio between the inside of the source region and the outside is smaller for 
long periods than for short periods, as expected. 

4. The gradient of intensity around the boundary of the source is steeper when the width of 
the source is wide. 
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5. Inside the source the ratio of non- 
cellular to cellular wave intensity is 
of the order of 10‘ 5 for r = 10 sec 
and 10" 3 for r = 100 sec. Although 
the none ellular intensity exceeds the 
cellular at large distances from the 
source, the contribution of the non- 
cellular wave to the observed in- 
tensity is negligible, because both 
waves attenuate in long distance 
propagation. 

6. According to the present calcula- 
tions, which are based on an 
isothermal atmosphere, the intensity 
of acoustic waves more than several 
hundred k m from the region of 
auroral activity is negligible. For 
the real atmosphere, therefore, the 
horizontal propagation of acoustic 
waves through the ducts around the 
mesopause and the stratopause is 
rather important in explaining the 
diurnal variation of the arrival di- 
rection of these waves during the 
period of high geomagnetic activity, 
as shown in Figure 1. 



7. 


It should be noted that auroral activ- 
ities are not necessarily within the 
auroral zone but, rather, extend 
toward lower latitudes (nearly to 


Figure 1 3— Intensities of cellular waves on the ground 
in the units of q 0 vs. the horizontal distance x from the 
center of the source, of which the width k Q is 100 km, 
f or r = 10 sec, 30 sec and 100 sec. Full lines and dashed 
lines stand for h = H and h = 1/2H, respectively. 


5 0°N), when these sonic waves are observed during the periods of magnetic disturbances 


(as can be seen from the auroral visoplot shown in Figure 2). 


To show the position and structure of sonic ducts as well as attenuation of the waves requires 
the actual atmospheric temperature distributions. This is discussed elsewhere (Reference 59). 

Finally, it should be noted that the energy flux of acoustic waves at the ground S in 
ergs/cm 2 -sec is 


S = Ec 


(78) 
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where 


p s and p s are the maximum amplitude of pressure change in dynes/cm 2 and the static density of 
air at sea level, and c is the sound velocity in cm/sec. 

Since p s is of the order of 1.25 x 10 " 3 gm/cm 3 , the energy flux corresponding to p s = 1 
dyne/cm 2 is approximately 1.4 x 10~ 2 erg/cm 2 -sec. To show the energy relation between the 
input power and the observed output intensity, as given in Figures 3, 5-10, and 13, the above re- 
lation must be used. 
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Appendix A 

Derivation of Differential Equation 37 


By differentiating Equations 30 and 31 with respect to t and using Equation 33 


d 2 u _ d , , > 1 ds 

- n \ c - -pis* • (Ai) 


d 2 w _ d , 2 » C dc 2_ l 1 as 

~ a^( c " [O'- ») e + srjx ~ p^ar 


(A2) 


When the time dependences of u, w, p, p, and s are proportional to a common factor e icrt , a 
straightforward calculation leads to Equations 34, 35, and 36 from Equations Al, A2, and 33, 
respectively. 

Differentiating Equations Al and A2 with respect to x and z, respectively, and adding them 
gives 


where 




1 ds 

P 0 C 2 <>* 



s 


(A3) 


dw du 
dx “ dz 


(A4) 


Similarly, by differentiating Equations Al and A2 with respect to z and x, respectively, and 
subtracting the former from the latter, 


d 2 £ _ dc 2 1 dy 1 /dc 2 \ ds 

= - pr + O'- 1) ej asr + pr + • (A5) 

The elimination of l from Equations A3 and A5 results in 



Since it is assumed that u, w, p , p, and s are proportional to e * CTt , the elimination of the time 
derivative from Equation A6 gives Equation 37. 
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Appendix B 

Solution of Equation 59 at z = 0 


Applying the Fourier transform to Equations 34 and 35 gives 


-a- 2 U = ike 2 X - ikgW - £7 S 

P o 


(Bl) 


-V 2 w - c 2 37 - ygX + ikgU - 37 


ik dS 
P 0 dz 


(B2) 


where U and w are the Fourier transforms of u and w . Solving these bilateral equations with re- 
spect to U and w yields 


( ct 4 _ k 2 g 2) u 


, dX f , , „ iker 2 „ ikg dS 

ikge 37 + ik( 7 g 2 - o' 2 C 2 ) X + 7^-(7y S + 37 


(B3) 


(cr«-k 2 g 2 )w - "c 2 cr 2 al + g(ycr 2 — k 2 c 2 ) X +7^7) s + Po CT (z) 3I 


k 2 


dS 


(B4) 


Since W(z = 0) = 0 and dS/dz z=0 = S( z =0)=0, from the above equations 


dX 

Hz' 


i-O 


- m(cr) X(O) 


(B5) 


where 


y k 2 

m ( cr > ° g iTT " 7T 


(B6) 


Similarly, the Fourier transform of Equation 33 is 


icrP(k, z) - gp 0 (z)W - c 2 P 0 (z)X(k, z) + S(k, z) 


(B7) 


where 


P(k, z) 


it 


e lkx p(x , z ) dx 
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Solution for X(z = 0), Direct Calculation 


From Equation 61 


dX 

dz 


C i - j F<y )e niZ dz' 

n 0 z 

+ n 2 e 2 

C 2 + ^7 J F(z ' ) e ”" 2 * ' dz' 



*0 


(B8) 


Using Equation B5, from Equation B8 we get 


n 2 - m(cr) 

" - _ — r* 

"l n 1 - m(cr) 2 


(B9) 


If we write 


A , = c 


r z 

i - ^7 J F(z') e'" 1 ' dz' 


and 


= C 2 + I" F (z') e " 2l dz ' 

** A 


then from Equations B8 and B9, by assuming h < 2H, 


Al e^ = 1 


Cj e' MZ for 0 < z < z 0 


o 4 - k 2 g 2 \ -( n ,^> 0 

C l +Q o 2(UC 2 CT 2 


P + k2 g _ k 2 

_ zl H -[(l/h)-(l/2H)]j 


+ Q 0 C 

0 2/r rij + £ 




(BIO) 


for z Q < z <oo 


and 


A 2 e^ z 


C 2 e^ z for 0 £ z < z 0 , 




[ c « - i 


k 2 g 2 1 -(n 2 +^)r 


2/i C 2 cr 2 n 2 + ^ 


— + k2 g - k 2 

eM . -0 - — e - [(!/»)-( 1/2H)]. 

0 2/x r> 2 + /3 


(B11) 


for z fl £ z < 00 , 
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where 


Qo 


i/2 sin^ 0 k (y-l)g 0 , o/h 
Vn k P s he 2 6 


(B12) 


and 


' (B13) 

The term including e represents the waves whose amplitudes increase with height. The cor- 
responding kinetic energies diverge to infinity provided that N 2 >M 2 . On the other hand, if N 2 <M 2 
this term represents inward- going waves. 

Because of the boundary conditions listed in the section ’’Mathematical Treatment,” the terms 
including e MZ must vanish. Therefore 


q 4 ~ k 2 g 2 


2m c 2 q 2 n 2 + 


-(n 2 +^)z 0 


(B14) 


Putting this constant into Equation 61, for z = 0, gives 


X(k, 0) 


n O'- ^qp sin\ 0 k g 4 - k 2 e 2 

\n h£T 4 c 2 Po e k (n, -m)(n 2 +/3) 


(B15) 


The continuity of W(k, z) and P(k, z) at z = z 0 is satisfied by the above choice of constants, 
i.e., it can be shown that 


limW^z 0 “ej - limW^z 0 +e'j 

limP(z 0 -e) = limPfzg+eM 
e-0 v ' e ' -*0 ' 


Solution for X(k,0] by means of the Laplace Transform 

By applying the Laplace transform with respect to z, 


X(k, s) = f e _sz X(k, z ) dz 


(B16) 


Equation 59 becomes 

, . . dX 

(s 2 - 2Ns + M 2 ) X - (s-2N)X(0) --37 


% 


- ( S+^ ) Z° 

s TJ j^s 2 - 2Ns - k 2 ( 1 - 



37 



By making use of the relation shown by Equation B5, this is 


(s-2N-m) „ Qq e 

X(k ’ s) " K - n \ II - n \ X(0) + 7 ? 

( s n ,)(s-n 2 ) 


(s+^)i 0 


T ** /* o r "i 

— 1- 7 n V 3 - 5 - . (B17) 

L ( s ~ n l) (s -n 2 )c 2 cr2 J 


where rij and n 2 are given by Equations 63 and 64; m = m(o-) and Q 0 are given by Equations B 6 
and B12. 


X(k, z) 


By making use of the inverse transform of Equation B16 the solution to Equation 59 is 

2 N - m - n, _ 2 N - m - n, \ 

— e" 2 '): 




2 m 


1 »i 

e 


2 m 


lx(0) + Q o e-f } '0(z-z o ) 


~ Qn 


4 ~ k 2 g 2 -0z o 
2/z C 2 cr 2 


n 2 (*"* 0 ^ n l(*”*o) 




|_ P + n 2 ^ + n i + (/3 + n,) (/3 + n 2 )J *(* Z °)’ < B18 ) 


Because of the boundary conditions at z = ®, the term proportional to e*** must vanish, i.e., from 
Equation B18 


2N - m - n 2 a 4 _ ^2 g 2 e ""(/ 3+n 2 ) I o 

2m X (°) ‘ % c 2 o ' 2 2 M(y 8 + n 2 ) = ° ' (B19) 

This gives the same solution as Equation B15. It should be noted that the above solution satisfies 
the conditions of continuity of w(x, z) and p(x, z) at z = z 0 , although u(x, z) and s(x, z) are not 
continuous at z = z 0 . 
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Appendix C 

Evaluation of p(x, 0) 


As can be seen from Equation B15 X(k. 0) is an even function of k. Therefore, the substitution 
of Equation B15 into Equation 67 gives 


x(*. 0) = — 


2 (y~ 1 >qo cr2 z 


,/2H 


hk 2 p c 4 

g ' S 


/•CO 


. cos xksin k 2 “ k 2 

„ J£ ,« 

0 k ( n i- m )( n 2 + ^) 


(Cl) 


where 


Since m is a real function of k for k >k c and complex for 0 <k <k c , where 


(C2) 


-4 


a 2 - a 2 


(C3) 


Equation Cl can be written 


x( x ’ 0) - (x, 0) + x nc (X, 0) 


(C4) 


where 


• ' ” hk.V.c- J, 


cos xk sin ^ n k k 2 - k 2 

o g 


( n i _m )( n 2 +P) 


dk , 


and 


X'nc ( X ’ IT 


2 (y“i)q 0 cr2 , 


hk 2 2 Ps c 4 


l /2 H 


f 


cos xksinX fl k k 2 “ k 2 
" 0 1 


( n l _m )( n 2 + ^) 


dk 


(C5) 


(C6) 


The integrand of y nc (x, 0) has two singular points &/ c and cr 2 /g, corresponding to the two 
types of free oscillations of the noneellular mode which were pointed out by Pekeris* and have 


•Pekeris, C. L., "The Propagation of a Pulse in the Atmosphere. II,” Phys. Rev. 73:145-154, January 15, 1948. 
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been discussed earlier. Because of a steep exponential term e M, ° , however, contributions of 
these singularities to the integral is not as important as the contribution from the narrow band 
near k >k c , where the exponential term is nearly unity. 

The integrand of x c has no singularities but oscillates by the term e IVZ °. The main contri- 
bution arises from a narrow domain near k < approx. k c . Equation 68 can be written 


p(x- 0) - P c (x. 0) + p nc (x, 0) 


(C7) 


where 


p c (*■ 0) = p r + ipi , 


(C8) 


Pc <*. 0) 


2 (y~ 

" hk^c 2 


V2H 


TV 11 ”' 

•>0 


cos xk sin k Q k 


k g 2 - k 2 


(n, -m) (n 2 + fi) 


dk 


= VP, 2 + P 4 2 . ( C 9) 


and 


, ~ 2<y-l)q 0 » , 0 /2H f k « 

P „ c(X ’ 0) ~ ▼ hk 2 c 2 6 J 


cos xk sin k 
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k 2 - k 2 
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( n l _m ) ( n 2 + / 8 ) 


dk 


(CIO) 
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